MULTI-HAMILTONIAN STRUCTURE FOR THE FINITE DEFOCUSING 

ABLOWITZ-LADIK EQUATION 



MICHAEL GEKHTMAN AND IRINA NENCIU 



Abstract. We study the Poisson structure associated to the defocusing Ablowitz-Ladik equa- 
tion from a functional-analytical point of view, by reexpressing the Poisson bracket in terms of 
the associated Carathcodory function. Using this expression, we are able to introduce a family of 
compatible Poisson brackets which form a multi-Hamiltonian structure for the Ablowitz-Ladik 
equation. Furthermore, we show using some of these new Poisson brackets that the Geronimus 
relations between orthogonal polynomials on the unit circle and those on the interval define an 
algebraic and symplectic mapping between the Ablowitz-Ladik and Toda hierarchies. 



It has been well-known since the work of Flaschka [Flalj , [Fla2 that the celebrated Toda lattice 
can be represented as an isospectral evolution equation for Jacobi matrices, by which we mean 
symmetric, tridiagonal matrices J, with positive off-diagonals. This not only allowed Flaschka 
and Henon to prove the complete integrability of the Toda equation, but also opened the door 
for the study of related Lie algebraic and geometric structures for generalized Toda flows with 
applications ranging from numerical analysis to quantum cohomology. We will not attempt to 
properly cite vast literature on the Toda lattice, instead referring the reader to the bibliography 
in |OPRS| or in the recent survey fBl oGek] . We will mention, however, two seminal works that 
introduced two aspects of the theory crucial to our exposition: Moser's paper [Mosj that explained 
the role played by spectral data in both explicitly solving and establishing complete integrability 
of the Toda lattice and Kostant's comprehensive treatment of the generalized Toda lattice as a 
restriction of the larger Hamiltonian system to a minimal irreducible coadjoint orbit of the Borel 
subgroup (the manifold of Jacobi matrices with fixed trace in the case of gl(n)) equipped with the 
Lie- Poisson structure. Combination of these two approaches later allowed Deift et al. to establish 
complete integrability of full Toda flows [DLNTlj . More recently, it was shown in |FayGek2] , that 
the spectral data can be used to define in a natural way a multi-Hamiltonian structure for a family 
of Toda- like systems associated with all minimal irreducible coadjoint orbits in gl(n). 
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In this paper we concentrate on the study of another integrable system, the defocusing Ablowitz- 
Ladik (AL) hierarchy, through its connection to a unitary analogue of Jacobi matrices, the so- 
called CMV matrices. The defocusing AL equation was defined in 1975-76 by Ablowitz and Ladik 
[AblLadlllAblLad2] as a space-discretization of the cubic nonlinear Schrodinger equation. It reads: 



(1.1) 



where /3={/3fe}cBisa sequence of complex numbers inside the unit disk and 
The analogy with the continuous NLS becomes transparent if we rewrite (jl.ljl aqj 

-0 k = Pk+i - Wk + Pk-i - \(3k\ 2 ((3 k +i + /3 fc -i). 
A simple change of variables, ak — e 2tt f3k for all k, transforms (jl.lj) into 
(1.2) - ia k = p\(a k+ i + a k -i), 



\Pk 



where pk = \/l — |afc| 2 . This is the equation we will refer to as Ablowitz-Ladik (or AL). In this 
paper we focus on the finite case, by which we mean the case in which a_i = —1 and, for some 
fixed n > 1, a„-i e S 1 = {z e C | \z\ = 1}. These are Dirichlet boundary conditions, and one can 
easily see from (jl.2|> that the evolution for ccq, . . . , a„_2 decouples from the evolution for the other 
a's. Consider the following Poisson bracket on the space of (ao, . . . , a,i_2, 0Cn-i) S D™ -1 x S 1 : 

df dg df dg 



(1.3) 



if, 9} 



N-2 
k=0 



da k da k da k da k 



where p k = yl — \oc k \ 2 and, for a complex variable a = u + iv, u, v e 
are defined as usual by 



I, the partial derivatives 



d_ 

da 



d_ 

du 



dv 



d_ 

da 



d_ 

du 



.d_ 

dv 



In this Poisson structure, the AL equation (|1.2p becomes completely integrable. Moreover, it can 
be re- written in the Lax form, with the Hamiltonian Re TrC, where the Lax operator C is the CMV 
matrix associated with the coefficients ao, . . . , a n -2, a n -i (for the background, see Subsection l2.1|) . 
In fact, one can define a whole hierarchy of evolution equations, that we will call the AL hierarchy, 
by considering the Hamiltonians given by the real and imaginary parts of K k = ^Tr(C fc ) for k > 1. 
In terms of Lax pairs, the hierarchy is given by the following evolutions equations (see jNenlj ): 



(1.4) 

and 

(1.5) 



{C, 2 ReOFT fc )} = [C, i{C k ) + + i((C k ) + )*} 



{C,2Im(K k )} = lC,(C k ) + -((C k ) + y 



for all k > 1, where for a matrix X, we have 



(X + )jk - 




if j < k; 
if j = k; 
if j > k. 



One of the central ingredients in the study of the AL hierarchy is, similarly to the Toda case, 
rewriting the Poisson bracket (jl.3|) as the restriction to the manifold of CMV matrices of the 
Gelfand-Dikij bracket on the associative algebra M„(C) of n x n matrices (for a short back- 
ground, see Subsection 12. 2\i . This was done independently by L. C. Li [EI], and R- Killip and 
I. Nenciu [KilNen2] . and allowed Killip and Nenciu to solve the system and find the long-time 
asymptotics of the as, and of certain associated spectral quantities. Inverse spectral problem for 
semi- infinite CMV matrices was recently utilized in this context in [Golj . Note also that an alter- 
native Lax representation was used to linearize finite and semi-infinite AL flows in [Gekj . while 
the approach based on continued fractions was suggested in [Com] . 



Here, and throughout the paper, / will denote the time derivative of the function /. 
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If dfi is the measure on the unit circle associated to the Verblunsky coefficients (ao, ■ • • j &n-2i &n-i) £ 
ID" -1 x S 1 , then it is known (see Subsection 12. ip that dfi — 2j=i Mj^zj . with fij E (0,1), 
\i\ ^ h /j, n = 1, and Zj — e l6j G S 1 for all 1 < j < n. The function 

n 

H Zj = dct(C) = (-l)™- 1 ^-! 

is a Casimir, and the manifold of CMV matrices with fixed determinant forms a symplectic leaf 
on which, for 1 < j, k < n — 10, we have 

(1.6) {0,- A} = ' to. I logK/Mn]} - fyfc, 
and 

(1.7) {log[^//x„],log[/x fe //i„]} = 2cot(^) +2cot(^^) +2cot(^ 1 ). 

These results are obtained by regarding the CMV matrices, or, equivalently, the associated spectral 
measures, as the central objects. In this paper, we adopt a slightly different point of view, and 
focus on the associated Caratheodory function, 

'C + z" 



F(z) 



C- 



li 



This is the analogue in the unitary case of the better known Weyl function (or m-function) 
associated to a Jacobi matrix, 

1 



m(A) = 



li 



In FayGck2 , Faybusovich and Gekhtman adopted this point of view for the Toda lattice and 
computed Poisson brackets induced by the Lie-Poisson structure for m(A) and ro(£) at any two 
distinct points A and £. The resulting Poisson structure on Weyl functions was then shown to be 
a part of a family of compatible Poisson brackets which constitutes a multi-Hamiltonian structure 
for the Toda lattice. In this paper, we follow the same road in the Ablowitz-Ladik case. 

The paper is organized as follows. In Section [5] we give some background information on 
the theory of orthogonal polynomials on the real line and unit circle, and on classical i?-matrices. 
Section[3]contains the first important results, Theorem[3]and Corollarv l3.3l which gives the formula 
for the Poisson bracket of the Caratheodory function at two distinct points in the complex plane. 
In particular, this represents a more direct proof of some of the results in [KilNen2] . Here we 
should also mention that in a recent paper [CanSim] , M. Cantero and B. Simon study Poisson 
brackets for orthogonal polynomials both on the real line and on the unit circle induced by standard 
Poisson structures for the Toda and Ablowitz-Ladik hierarchies, respectively. An essential part of 
their analysis are the formulae for the Poisson brackets of Weyl and Caratheodory functions, for 
which they give new proofs by induction, using purely orthogonal polynomial methods. Another 
related recent paper is [Tsij , which also introduces a family of Poisson brackets compatible with 
the Sklyanin bracket associated with the standard 2x2 rational solution of the Classical Yang- 
Baxter equation. These Poisson brackets are defined on monodromy matrices associated with 2x2 
spectral parameter depending Lax representation for a family of integrable systems that includes 
both open and periodic Toda lattices. 

The formula (|3.12[) for the Poisson bracket mentioned above allows us to extend the Poisson 
structure to the space of finite, but unnormalized measures on the circle, and obtain the canonical 
coordinates for both the extended and the usual Poisson structures; this is achieved in Section 2J 
Finally, we define the family of compatible Poisson structures in Section and show its connection 
to the defocusing Ablowitz-Ladik equation. In the last section, Section [51 we go back to the Toda 
lattice and show that, loosely speaking, 'half of the Ablowitz-Ladik hierarchy (also known as 
the Schur flows) is mapped symplectically onto the Toda lattice hierarchy via the well-known 
Geronimus relations. 



2 Since the ZjS are always distinct, any choice of labeling for z\, ...,Zn is locally well-defined and leads to these 
formulae. 
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2. Background 



2.1. Orthogonal polynomials. As CMV matrices arose in the study of orthogonal polynomials, 
it is natural that we begin there. We will first describe the relation of orthogonal polynomials to 
Jacobi matrices and then explain the connection to CMV matrices. 

Given a probability measure dv supported on a finite subset of R, say of cardinality n, we can 



apply the Gram-Schmidt procedure to {1, x, 
L 2 {dv) consisting of polynomials, {pj{x) 
this basis, the linear transformation f{x) 



,x" } and so obtain an orthonormal basis for 
j = 0, . . . , n — 1}, with positive leading coefficient. In 
-> xf(x) is represented by a Jacobi matrix, 



(2.1) 



with a,j > 0, bj e M 
three-term recurrence: 



J 



a i 



bn 



An equivalent statement is that the orthonormal polynomials obey a 



xpj(x) = a,jPj+i(x) + bjPj(x) + a,j-ipj-i(x) 



where a_i = and p n = 0. A third equivalent statement is the following: A is an eigenvalue of J if 
and only if A 6 supp(d^); moreover, the corresponding eigenvector is [po(A),pi(A), . . . ,p„_i(A)] T . 

We have just shown how measures on M lead to Jacobi matrices; in fact, there is a one-to-one 
correspondence between them. Given a Jacobi matrix, J, let dv be the spectral measure associated 
to J and the vector e\ = [1, 0, . . . , 0] T . Then J represents x i— > xf(x) in the basis of orthonormal 
polynomials associated to dv. 

Before explaining the origin of CMV matrices, it is necessary to delve a little into the theory 
of orthogonal polynomials on the unit circle. For a more complete description of what follows, 
the reader should turn to Sirnl]. Given a finitely-supported probability measure dp, on S , the 
unit circle in C, we can construct an orthonormal system of polynomials, <f>k, by applying the 
Gram-Schmidt procedure to {1, z, z 2 , . . . }. These obey a recurrence relation; however, to simplify 
the formulae, we will present the relation for the monic orthogonal polynomials $fc(z): 



(2.2) 



4> 



fc+i 



Here ctk are recurrence coefficients, which are called Verblunsky coefficients, and denotes the 
reversed polynomial: 

k k 
(2.3) * fc (z) =^>z< => *l(z) =Y,Ck-iz l = z k <S> k {l/z). 



1=0 



1=0 



When d/i is supported at exactly n points, a k 6 B for < k < n — 2 while a„_i is a unimodular 
complex number. (Incidentally, if dfi has infinite support, then there are infinitely many Verblun- 
sky coefficients and all lie inside the unit disk.) The Verblunsky coefficients completely describe 
the measure dfi: 

Theorem 1 (Verblunsky). There is a 1-to-l correspondence between probability measures on 
the unit circle supported at n points and Verblunsky coefficients (a , . . . ,a„_i) with a k £ B for 
< k < n — 2 and a n _i G S 1 . 

^From the discussion of Jacobi matrices, it would be natural to consider a matrix representation 
of f(z) i — > zf(z) in L 2 (dp). Cantero, Moral, and Velazquez had the simple and ingenious idea 



to define a basis in L {dp) by applying the Gram-Schmidt procedure to {l,z, 



}■ 



The resulting functions, Xk{z) (0 < k < n 
polynomials: 



1), are easily expressed in terms of the orthonormal 



(2.4) 



XfcO) 



-(fc-l)/2 



(f>k{z) 



k even 
k odd. 
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In the orthonormal basis {xk( z )} of L 2 {dp), the operator f(z) i— > zf(z) is represented by the 
CMV matrix associated to the Verblunsky coefficients of the measure d/i: Given the Verblunsky 
coefficients cxq, . . . , a n -i in D and o„_i 6 S 1 associated to the measure d/i, let pk = \/l — |a/c| 2 , 
and define 2x2 matrices 



for < k < n — 2, while 5_i = [1] and 3 n _i = [a n _i] are I x I matrices. From these, form the 
n x n block-diagonal matrices 

C = diag(S ,S 2 ,S 4 , . . .) and M = dia I g(2_i,H 1 , 3 3 , . . . ). 

The CMV matrix associated to the coefficients ao, . . . , a n -i is C = £AL 

The measure dfi can be reconstructed from C in a manner analogous to the Jacobi case: 

Theorem 2. Let dfi be the spectral measure associated to a CMV matrix, C, and the vector e\. 
Then C is the CMV matrix associated to the measure dfi. 

Proofs of these Theorems can be found in [CanMor Vellj or [Simlj . As explained in the In- 
troduction, throughout the paper we will always use implicitly the bijection between measures 
dfi = J^^Mi) CMV matrices, and the coordinates given by the Zj's and fij's. A very impor- 
tant notion, that will be heavily used in this paper, is the Caratheodory function associated to a 
probability measure fi on the unit circle S 1 ; it is given by (see [SimH Section 1.3]) 

(2.5) F(z)= [ f£±±dnV). 



In terms of the other coordinates, F is given by 



C Z ' 11 j—i Z 3 



where C is the CMV matrix associated to the measure fi. Moreover, the Caratheodory function is 
related to the Schur function f by 

1 ' W l-zf(z) JK ' z l + F(z) 

While the Caratheodory function plays a very important role throughout the theory of orthogonal 
polynomials on the unit circle, it is very simple to see from the formulae above that, in the case 
of finite measures, F encodes exactly the same information as the measure /i. This is exactly the 
reason why we can use the Caratheodory function in a functional analytic approach to describe 
all the Poisson structures that we will introduce. 



2.2. Integrable systems and classical R-matrices. The manifold of Jacobi matrices with fixed 
trace forms a co- adjoint orbit of the group of invertible upper triangular matrices, if one views 
the space of symmetric matrices as a dual space to the algebra of upper triangular matrices. The 
Lie- Poisson structure on this dual thus induces a symplectic structure on the Jacobi orbit. More 
generally, the space of 3-diagonal (not necessarily symmetric) matrices form a Poisson submanifold 
in gl(n; R) with respect to the Lie- Poisson bracket associated to a particular Lie algebra structure 
on the n x n matrices (although not the one defined via the usual matrix commutator). These 
matters are described in detail in [Dei( IOPRS|, IPerj . for example. In contrast, CMV matrices are 
elements of the unitary group and hence the natural backdrop for CMV is that of Poisson-Lie 
groups or, more specifically, the group Gl(n,C) equipped with the Sklyanin bracket (see [OPRS. 
§2.12]). However, we choose to give a presentation in which the associative algebra of n x n 
matrices takes center stage; an analogous construction for KdV using the algebra of pseudo- 
differential operators was given by Gelfand and Dikij |GelDik| . This approach is described in 
Section 2.12.6 of |OPRS| . (Note that here we are referring to the second symplectic structure 
associated with KdV, which was originally proposed by Adler |Adl|. §4].) 
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Let g denote the (associative) algebra of n x n complex matrices. The algebra structure gives 
rise to a natural Lie algebra structure: 

[B, C] = BC- CB. 

As a vector space, g = i © o, where 

a={A: A=-A*}, 

is the space of skew-Hermitian matrices, which is the Lie algebra of the group U(n) ofnxn unitary 
matrices, and 

[={4eg: Lij = for i > j and L ifi G R} 

is the space of upper triangular matrices with real diagonal entries (the Lie algebra of the group 
L(n) ofnxn lower triangular matrices with positive diagonal entries). We will write 7r and 7T( 
for the natural projections onto these summands. This vector-space splitting of g gives rise to a 
second Lie algebra structure. First we define R : g — > g by either 

R{X) =tti(X) -ir a (X), for all X eg, or 

' " " f R(L + A) = L - A for all A e a and L e I 

The second Lie bracket can then be written as either 

[X,Y] R = ±[R(X),Y} + ±[X,R(Y)] VX,Y£g, or 
"" ; [L + A,L' + A'] n = [L, L'\ — [A, A'] V L, U e I, and A, A' e a. 

The second definition also makes it transparent that the bracket [ , ] R obeys the Jacobi identity. 
But the main property that we are interested in is that R obeys the modified classical Yang-Baxter 
equation: 

[R(X),R(Y)]-R([R(X),Y] + [X,R(Y)})=-[X,Y] 

This allows us to define Poisson brackets as follows: We can identify the dual space g* with g 
using the pairing 

(2.10) (X,Y) =ImTr(XF). 
The form ( , } is is non-dcgcncratc symmetric and invariant: 

(2.11) (X,[Z,Y\) =ImTr(XZY - XYZ) = ([X,Z],Y) 
or, equivalently, 

(2.12) (BXB~ 1 ,BYB~ 1 ) = (X,Y), for any B e GL(n, C). 
Given a smooth function ip : g — > R and B G g, define Vip : g — > g by 

d 



v(B + tC) = {V<p\ B ,C). 

t=o 

Equivalently, if we write bk,i — Uk,i + ivk,i for the matrix entries of B, then 

(2.14) mkil= **L +i »(L. 

We can now define the desired Poisson brackets on g. Let R : g — > g and (• , •) be as above. 
Given ip t , ip 2 '■ Q — > let Vj = S/ifj for j = 1, 2. Then both 

(2-15) Wi,V2} LP \ x = l([V 1 ,V 2 } R ,X) 

and 

(2.16) Wi,^2} GD \ x = \{R{^iX), V 2 X) - ±(fl(XVi), XV 2 ) 

define Poisson structures on g, the first one known as the Lie-Poisson (LP) bracket, and the 
second known as the Gelfand-Dikij (GD) bracket. Then under the embedding J i— > iJ, the 
manifolds of Jacobi matrices with fixed trace are symplectic leaves in (g, {■, -} LP ); this is just 
the usual construction in Gl(n, R) in disguise. Similarly, the manifold of CMV matrices with 
fixed determinant forms a symplectic leaf in the Poisson manifold (g, {-, -} GD ). Furthermore, the 
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restriction of the GD Poisson structure to the manifold of CMV matrices coincides with the AL 
Poisson structure (|1.3p (see Ki lNen2| ). 



3. Poisson brackets for Caratheodory and Schur functions 

Let Vi and V2 be the gradients at a point X of two functions (pi and ip 2 , as defined in 
Subsection 12.21 Then a simple calculation shows that the two brackets defined above can be 
written in the following slightly modified form: 

(3.1) {<p x , <p 2 } LP = \([X, Vj,ii(2V 2 )> - \({X, V 2 ],iJ(2Vi)> 
and 

(3.2) {<pi,<f»} aD = i<[X, W 1 ],R(XV 2 + V 2 X)) - \([X, V 2 ],R(XV 1 + ViX)>. 

Before going any further, we give the following definition. Consider k > and two smooth 
functions ipi and tp 2 , with gradients Vi and V 2 respectively. Then we set 

(3.3) 4{v?i, (p 2 } (fc) = (x, [Vi, i?(X fc V 2 + V 2 X k )] + [i?(X fc Vi + V^), V 2 ] ). 

One must note that this expression is a Poisson bracket only for k = 0, in which case it is exactly 
the Lie-Poisson bracket, and for k = 1. In this latter case a simple calculation using the properties 
of R and (|3.2|) shows that {•, is actually the Gelfand-Dikij bracket defined in (|2.16|) . But even 
though (|3.3p is not in general a Poisson bracket, in what follows we will work with the general 
expression as it makes it easier to emphasize the relevant steps of our calculations. 
We now focus on functions given by 

<p(X) = ImTr(P$(X)P) = Im($(X)u), 

where $ is a smooth function on M n (C), and P denotes the orthogonal projection on the vector 
ei = (1,0,...,0) T , 

P = e\e 1 . 

Note that any (p defined this way is invariant under the conjugation by invertible matrices with 
off-diagonal entries in the first column and row all equal to zero: 

tpiCXC- 1 ) = ip(X) 

for 

c =(o I 

where c\ 7^ and C 2 is an invertible (n — 1) x (n — 1) matrix. 

Lemma 3.1. Let ip(X) = Im($(X)n) as above. Then for any matrix B € M n (C) we have 

(3.4) (V(p\ x ,[X,B]) = (Vcp\ xl {X,BP + PB}) = {lV<p\ x ,X},BP + PB) 
where, as above, P — efei is the orthogonal projection on the vector ei = (1, 0, . . . , 0) T . 
Proof. Let C = exp(iA), where A is any matrix of the same form as C . Then 



dt 



v(CXC- 1 ) = {Vp\ x ,{A,X}) 



which implies that, for any B an expression (Vip\ x , [X, B]) does not depend on the (n— 1) X (n— 1) 
submatrix of £? obtained by deleting the first row and column. □ 
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Using the properties of the P-matrix and Lemma 13. II we see that: 

4{v>i, <M (fc) = ( [X, Vi] , R(X k V 2 + v 2 x fe )^ - / [X, v 2 ] , R{x k y x + \7 1 X k ) 
Vi, [P(X fe V 2 + V 2 X fc ),X]) - (v 2 , [i?(X fe V! + V!X fc ),X] 
Vi, [PR{X k V 2 + V 2 X k ) + R{X k \7 2 + V 2 X k )P, X] } 
- /v 2) [Pi?(X fc Vi + ViX fe ) + R{X k Vi + ViX k )P, X] 
= ( [X, Vi] , PR(X k V 2 + V 2 X k ) + R(X k V 2 + V 2 X fe )p) 
[X, V 2 ] , PR{X k Vi + ViX fc ) + P(X fe Vi + ViX k )P 



While this expression appears to be more complicated than the one we started from, it will shortly 
be shown that it is exactly what we need in order to continue our calculation. 

Using Lemma ETTl we can compute the Poisson bracket for functions tp as above. Note that this 
calculation is general and extremely robust, and it applies to a large range of expressions involving 
P-matrices. The first such calculation that we are aware of was done for the finite Toda lattice 
by Faybusovich and Gekhtman FayGek2 . Proposition 13. 21 gives a short proof of the next step in 
the calculation, while Theorem [3] takes it to its conclusion. 

Proposition 3.2. Let ifj(X) = Im($j(X)n) = ImTr(<Fj(X)P) for j = 1, 2, where 3>j are smooth 
functions, and set Vj = Vtpj(X). Then the values of the linear and quadratic brackets of the tpj 's 
are given by: 

(3.5) 2{ Vl , V2 } LP \ x = Im(2(XV 1 V;) 11 - 2(XV 2 VJ)n + [X, [Vi, V 2 ]] n 
and 

2{V1. V2\ GD | x = I™h(XV 1 V* 2 X*) n - 2(V 1 XX*V* 2 )n 

(3.6) V 

+ [XV 1 ,XV 2 }n-[ViX,V 2 X] 11 



Proof. First note that for any matrix A we have 

P R( A) = 2Re( Au)P - PA and R(A)P = AP + 2A*P - 2Re(A n )P. 

Therefore 

PR(A) + R{A)P = [A, P] + 2A*P 
and so for any fc > the expression 4{t^i, </? 2 }^ equals 

( [X, Vi] , 2(X k V 2 + V 2 X k )*p) - ( [X, V a ] , 2(X k V 1 + ViX fc )*P^ 

( 3 - 7 ) / v / 

[X, Vi] , [X fc V 2 + V 2 X fe ,P] ) - ( [X, V 2 ] , [X fe Vi + \7,X k , P] 



Taking (|3~7|) for fc = and 1 yields the result. □ 

Note that the P-matrix does not appear any more in (|3.5[) or (|3.6[) . This allows us to prove 
our first main result: compute the respective brackets for the Weyl and Caratheodory functions. 
More precisely, let X € M n (C) be a matrix, and z, A € C \ spec(X) complex numbers. For k > 0, 
consider 

(3.8) mW(A) = (X fc (A - X)- 1 ) n = Tr(X fc (A - X)~ l P) 
and 

(3.9) P w (z) = (X k (X + z)(X - z)- 1 ) u = Tr(X k (X + z){X - z)- x P). 

We think of these functions as being defined on the space of matrices, and depending on a complex 
parameter, A or z respectively. Note that, for fc = 0, we recover the Weyl and Caratheodory 
functions, respectively. We then have the following result: 



MULTI-HAMILTONIAN STRUCTURE FOR AL 



9 



Theorem 3. Let X £ M„(C) be a fixed matrix, and z,w, A, £ 6 C\spec(X) 6e complex parameters, 
with z ^ w and A ^ £. Then, for k > 0, 

m(A) — m(£) 



(3.10) {m(A),m(0} (fe) = (m (fe) (A) - m (fe) (£)) 
and 

(3.11) {F(z),FH} (fc+1) = i(F<*>(z) - fW(w)) 



m(A)m(£) 



F(z)F(w) - 1 - ^-±^(F(z) - F(iu)) 
z — w 



Corollary 3.3. LeiC be a finite CMV matrix, and F and f the Caratheodory and Schur functions 
associated to the spectral measure of C and e\ = [1, 0, . . . , 0] T . Then, for two distinct points 
z, w £ C \ spec(C), the GD brackets of these functions at C are given by: 

{F(z),F(w)} GD = i(F(z) - F(w))(F(z)F(w) - l) 

( 3 ' 12 ' ) .Z + W , ,2 

— i [b [z) — b (w)) 



an 



d 



(3.13) {/(*), /H} GD = -2i f{z) z _ f w {w) {zf(z) wf(w)). 

Proof. The first relation is just a special case of p. lip for k = and X = C. The bracket (|3.13[) 
follows from (12.71) and the observation that 



1 d 


n-F( z )\ 




a-F(w)\ 


zw dF(z) 


\l + F(z)J 


dF(w) 


\l + F(w)J 



{f(z),f(w)} GD = —-^ ( T - 7 A-( ) ( - { ^ ) {F(z),FH} GD 

1 



(l + ^(z)) 2 (l + ^(w)) 2 

Now use (|3.12p and the expression (|2.7p of F in terms of /. □ 

Proof of Theorem [3J The proofs of both relations follow the exact same ideas, but since in this 
paper we focus on the Ablowitz-Ladik system, and hence the Gelfand-Dikij bracket, we will only 
give the proof of p. lip . We approach this by first computing the fc-brackets of Im F and Re F = 
Im(iF) at two different points z and w. 

Let R be the R- matrix defined in Section [5J We start by working with the general expression 
for the fc-bracket and for any two functions (fx and ip 2 as above. We know that 



4{^i, <^ 2 } (fc) (X) = { [X, Vi] , PR{X k V 2 + V 2 X fe ) + R(X k V 2 + V 2 X k )P 
- ( [X, V 2 ] , PR{X k V l + ViAT fc ) + R{X k V l + ViAT fc )P 



Just as in the proof of Proposition 13. 21 the observation that allows us to continue is that, for any 
matrix A: 

PR(A)=2Re(A u )P-PA and R(A)P = AP + 2A*P - 2 Re(A u )P. 

Therefore 

PR{A) + R{A)P = [A, P] + 2A*P 

and so 



(3.14) 



4{^ 1; V2 } (k) {X) = ( [X, Vi] , 2(X k V 2 + \7 2 X k )*P 

- ( [x, v 2 ] , 2(x fc Vi + y x x k )*p 

'[X,Vi], [X k V 2 +V 2 X k ,P] 
'[X,V 2 ],[X k V 1 + V 1 X k ,P] 
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We wish to apply this formula to the real and imaginary parts of the Caratheodory function F 

'X + z" 



11 



Indeed, let u(z) = ReF(z) = ImiF(z) and v(z) — ImF(z). Note that, for a fixed parameter 
z, these functions are well-defined in a neighborhood of X, and they are of the type we have 
considered above. If we denote V z = Vv(z), then we get 

Vu(z) = iW z = -2iz(X - z)- 1 P{X - z)- 1 

and 

Vv(z) = V Z = -2z(X - z^PiX - z)- 1 . 
Plugging these into the formula for the k th i?-bracket, we get 

{u(z),u(w)}W = Im(T 1 (fc)-T 2 (fc)-T 3 (fc)+T 4 (fc)) 
{v(z),v{w)}W = Im(Ti(fc)-T 2 (fe)+T3(A;)-T4(fc)) 
{u(z),v(w)}W = Im(^^ 1 (fc)+^^ 2 (fc)+^^ 3 (fc)-^^4(fc)) 

= Re(Ti(fc) + T 2 (k) + T 3 (k) - Ti(k)) 
{v(z),u(w)} (k) = Im(-*ri(fe) - iT 2 (k) + iT 3 (k) - iT 4 (k)) 
= Re(-T 1 (fc)-T 2 (fc)+T 3 (fc)-T 4 (fc)), 



where 



Ti(k) 
T 2 (k) 
T 3 (k) 
Ti{k) 



iTr(2[X,V 2 ] (X k V w + V w X k )*P 
\Tr(2[X, V w ] (X k X7 z + V z X k )*P) 
^Tr([X,V z ] • [X k V w + V w X k ,P]) 
±Tr([X,V w ] ■ [X k V z + V z X k ,P]). 



Therefore 



{F{z),F{w)} {k) = {u{z),u{w)} {k) - {v(z),v(w)}W 

+ t({u(z),v(w)}^ +{v{z),u(w)}^) 
= -2Im(T 3 (fc) -T 4 (fc)) + i-2Rc(T 3 (fc) - T 4 (fc)) 
= 2i(T 3 (k)-T 4 (k)) 

Let us note in passing that in the k = 1 (Gelfand-Dikij) case, we have that T\ = T 3 and T 2 = T 4 . 

In order to compute T 3 — T 4 , we also need the following simple observation: For any three 
matrices A, B and C, 

(3.15) Tt(APBPCP) = Tr(PAPBPCP) = A n B n C n - 

After simplifying and grouping terms together, we get that 

4T 3 (fc) - 4T 4 (fc) = I(fc) + II(fc) + Ill(fc), 

where 

I(fc) = 8zw((X z)- 1 X k (X W ;)- 1 ) n .[(X(X - - 

-(X(X-w)- 1 ) 11 ((X-z)-i) 11 _ , 

II(fc) = 8zw((X z)-'X{X W )- 1 ) 11 .[(X fc (AT - z)- 1 ) n ((X - «,)"% 
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and 



Ill(fc) = 8zw({X - z)~ 1 (X - w)- 1 )^- (X(X - z)- 1 ) 11 (X k (X - w)- 1 ) u 

-(X(X-w)-') u (X k (X-z)-') 11 



j ez. 



Recall that, for j > 0, wc have defined the function 

(3.16) F U) {z)= ^^±l 

Then we have that, for j G Z, 

(3.17) (X^(X z)- 1 ) n = \(F^(z) + F«(0)) - Yz^ FU+1){z) ~ fU+1) ^)> 



((X-z)- 1 X^(X-w)- l ) n 



1 



((X-z)-\X-w)- 1 ) 11 



1 



2(z - io) 



(F^(z)-F^(w)), 



1 



(F(z) - 1) - -(F(w) - 1) 



(3.18) 
and 

(3-19) ,,, 

11 2(z — w) \z w 

We will use the appropriate formula of ()3.17[1 — (|3.19|) in order to express our quantities only in 
terms of F = F^ and F<~ k \ Thus we get 

(z-w)- I(k + 1) = (F^(z) - F^(w)) [z(F(z) + l)(F{w) - 1) 

-w(F(w) + l)(F(z)-l)], 



(3.20) 



(3.21) 

and 
(3.22) 



(«-«>)• H(A + 1) = (F(z) - F(w)) [z(F^(z) + F^(0))(F(w) - 1) 

- w(F^ (w) + F^(0))(F(z) - 1)] , 

(z-w)- lll(k + 1) = (w(F(z) - 1) - z(F(w) - 1)) x 



[(F(z) + l)(fWW+fW(0)) - (f (w) + l)(FW(z) + f(*'(0))]. 
A straightforward calculation shows that II + III = I, and hence we find that 

(3.23) {F(z),F(w)}( k+ V = i(f (*>(*) - F^H) [f(z)F( w ) - 1 - - F(w)) 

as claimed. 



z — w 



□ 



We wish to use (13. 12[) to find the bracket of the eigenvalues and masses of the spectral measure, 
fi = Y^j=\ fizjUj, z j = e lBj and Yl Mj = 1, of a CMV matrix C. Note that, since for the next couple 
of sections we only work with the GD bracket, we will not specify it in order to simplify notation. 
We hope that this will not cause any confusion. 

Fix a CMV matrix C, with d[i = X^-=i ^z^j, Zj = e l6j and = 1, the associated spectral 

measure. Let us expand the bracket {F(z), F(w)}\ c in terms of the zjs and /ijS: 



" r 

{F(z),F(w)}= 

j,k=i 



Z + Zj 



Z — Zj 



W + Zfe 



j W - Zk 

Z + Zj W + Zfe 



{/J-j,vk} 



n 

V 

Ew + Zk d ( z + Zj 
H ~r~ 



{Zj,Vk} 

•'k «*j \ « *j / 

d / z + Zj \ <i / iu + Zfe 
dz 7 - \ z — Zj ; y dzfe V w — Zfe 



{Zj,Zfe.}. 
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Since 

d (z + C\ _ 2z 

we obtain 

rev \ T?f M (Z + + z fc) r , 

{F{z),F(w)}= > ^ 

^i (*-*i)(«>-«k) 

A 2z{w + z k )nj 

ki ^-^-^) 2 ' {Z3,Zkl 

Let 1 < s, t < n. If s ^ i, then we choose T s and T t to be small, positively oriented contours 
around z s and z t , respectively; we require that they do not intersect, nor surround more that one 
eigenvalue. In the case s = t, we choose two small, positively oriented contours T s and r' s around 
z s so that the contour T s is completely contained in the interior of T' s . By the residue formula the 
previous expansion implies that 

(3.24) aAi / [ {F(z),F(w)}dzdvj = 4z,Zt{n a) in], 

(3.25) ,J~-\o f f {z - z s ){F(z),F(w)}dzdw = Az s ztn s {z s ,fi t } 7 
{2m) 2 J Tt J Fs 

and 



(3.26) (2mY ' ' ( z ~ z ^)(w - z t ){F(z), F(w)} dzdw = 4z s z t ^ s ^ t {z s , z t }. 



Note that, if s = t, then we set T t — T' s ; in other words, we first integrate over the smaller of the 
two contours around z s . 

Theorem 4. With the definitions from the previous sections we have that, in the GD Poisson 
structure, 

(3.27) {z s ,z t } = 0, 

(3.28) {z s , = 2iz s pL t {8 st - Ms), 
and 



(3.29) {^siVt} = 2i\i s \i t 



Zk + z s sr^ z t + z k z s + z t 



7 - L -s Zk ~ Zs k~^t Zt ~ Zk Z s — Z t 



Proof. We prove the result by finding the residues generated by the right-hand side of 

{F(z),F(w)} =i(F(z) - F(wj) (F(z)F(w) - l) 

z — w 

in the integrals from . (|3~23]) . and (pOS)) . 

We begin with (|3.26[) . Note that only the quadratic poles in both z — Zj and w~z k play any role. 
But for s ^ t there is no term on the right-hand side of (|3.30[) which contains the denominator 
(z — z s ) 2 (w — z t ) 2 , and hence the double integral over T t and T s is identically 0. This proves (|3.27[) . 

Now we turn to (|3.25p . Assume first that s =fc t. By the same reasoning than above, the 
only terms in (|3.30[) which contribute to the integral are the ones containing the denominator 
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(z — z s ) 2 (w — z t ). In this case, that translates into 



1 



(2wi) 2 
1 

(2^)2 



r t Jr 



(z - z s )iF(z) 2 F{w) dzdw 



r t Jr 



(z + z s ) 2 2 w + z t 
l\z — z s j- r^fJ-s A 4 * dzdw 



(z - z s y z t - w 



{-i)Azifxi2ztfit 



or, equivalently, 
(3.31) 



{z s ,fi t } = -2iz s \x s \x t 



for s 7^ t. The case s = t is slightly more complicated because the factor (z + w)/(z — w) in the 
second term of the right-hand side of (|3.30|) plays a role. Indeed, 



1 



(z - z s ){F(z), F(w)} dz = i4z 2 sf , 2 F(w) - i^^Az 2 ^. 



Recall that T s is contained in the interior of F' s , and hence the function z — > is analytic on 
an open neighborhood of the interior of T s as long as w is on I^. Now integrate on I^: 

4:Zs/J, s {z s ,fi s } = -8iZs(x 2 z s n s - i(-2z s )4ZsfJ, 2 = Az 2 ^ s (-2iz s ^ 2 s + 2iz s fi s ). 

Combining this with (13.31[) gives (|3.28p . 

Finally, we turn to (|3.24[) . Here the formulae are more involved and all the terms play a role. 



^- jf {F(z), F(w)} dz = -2iz si i s (f(w) £ ■ 



- 2iz s [i s F(w) I ^2 



k^s 



Z s + Zfc 
Zk' 

z s + Zk 
Zs %k 



[M k - F(w) 



4iz s fi s 



w 



E 



k^s 



Zs + Z k 



fik - F(w) 



Integrating this on T t for t ^ s yields 

4z s Zt{fj, S) n t } 



siz s z t ^ z ^ t y 



fc#s 



8iz s z t ^ s fi t ^2 



Zs + Zk 
/ 

Zk 

Zt + z k 



8iz s z t ^ s ^t 



kjtt 
Z s + Z t 
Zs ~ Z t ' 



z t - z k 



or, after simplifications, 

{Ms,Mt} = 2[i s Ht 

as claimed. 



fc^s 



Zk 



Zk 



-Ilk 



E 

fc/t 



. Zi + Z k Z s + Z t 

I Hk + 1 



Zt - Zk 



Zs — Zt 



□ 



Remark 3.4. As an immediate consequence of Theorem|4]we can recover the analogous formulae 
for the Ablowitz-Ladik system obtained by Killip and Nenciu, |KilNen2] . Indeed, let us rewrite 
the brackets in terms of 6 S (where z s = e l6a ) and log[/it]. We easily obtain from (|3.27| and (|3.28l) 
that 

{0 s ,9 t } = 

and 

{9 S) log[fit\} = 25 st - 2[i s . 
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The first formula is a direct consequence of the Lax pairs in |KiiNen2, Proposition 4.5], while the 
second one coincides with equation (72) from the same paper. Consider (|3.29p and use the fact 
that 

l— : = COt I 

Then we get that 



ri r l r n X - ^ / @k — @s \ ^— \ / Of — Ok \ ( 9 S — 9t\ 

{log[/^],log[/x t ]j = }^ cot { 2 — J /ifc + l^ C0t { — 2 — + cot { — 2 — ) 



k^s k^t 



E 



< <>t( ® k - j col ! — j + col 



Mfc- 



k^sj L 

In the last identity we use the fact that 

1 - n s - fM = ^ t ik - 

k^s,t 

Note that we recovered Proposition 8.5 from [KilNen2] . 

While our proof of (|3.27| and (|3.28p is not necessarily shorter than the one from [KilNen2] . 
the case of (|3.29j) is completely different: not only is this derivation much simpler and shorter, it 
also only uses the i?-matrix formulation of the Gelfand-Dikij bracket. The approach of Killip and 
Nenciu exploits the asymptotics of the spectral parameters, and the expression of the bracket in 
terms of Verblunsky coefficients. 

4. The extended bracket for unnormalized measures 

We will now extend the bracket from probability measures /i on the unit circle 5* 1 to general, 
finite measures (~i on S 1 . If c = ^(S 1 ) is the total weight, and F the Caratheodory function 
associated to the normalized measure, note that c and F fully characterize the unnormalized 
measure. In other words, in order to define the extended bracket it is sufficient to give its value 
for these two quantities: 

(4.1) {c 1 F(z)} = ic(F(z) 2 -l) 1 

and then extend it to all smooth functions using the Leibnitz rule and bi-linearity0 While this 
describes the bracket uniquely, it is not clear that it obeys the Jacobi identity. In order to show 
this, we will rewrite the bracket in a different set of coordinates. Indeed, from (|4.1[) we get, by the 
same methods in the proof of Theorem [U that 

{c, z s } = -HcZsHs and {c, fx s } = -2icfj, s ^ s + Zj , 

where for a finite measure jl we always denote by fi the associated probability measure, /i = /Et/|/i| 
and, as in the previous sections, [i = Y^j—\ Mi'^zji z j — e *• Further set flj = c/j,j to be the 
corresponding weights for the unnormalized measure /i0 Note that the space of measures jl is 2n- 
dimensional and parameterized by 9\, . . . , 8 n , . . . , fl n . In these variables, a direct calculation 
using (|3.27p - (|3.29p and (|4.ip shows that for all 1 < s, t < n the bracket is given by 

(4.2) {8 Sl 6 t } = 0, {6 s ,fi t } = 2(M t 5 st , 
and, for s =/= t, 

(4.3) {As, At} = 2ip>sfk— — : — = 2/t s /t t cot 

z s — Zt 

So now we can easily prove that 



■^Note that here we use the same notation for the Poisson bracket on the space of finite measures as for the one 
on the space of probability measures. But since the later is just the restriction of the former, we trust that this will 
not create any confusion. 

^Recall that any symbol with a tilde refers to the unnormalized measures. 
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Proposition 4.1. The bracket given by (|4.2|) and (14. 3|) on the space of finite measures supported 
at n points on the unit circle obeys the Jacobi identity, and hence it defines a Poisson bracket. 
Furthermore, this bracket is nondegenerate. 

Proof. Checking the Jacobi identity is actually a simple calculation if one uses the following (slight) 
variations of (|4.2[) and (|4.3p obtained if we replace the jljS by | logfAj]: 

{0 s ,d t } = 0, {6 s ,llog\fi t ]} = S st , 

and 

{ilog[£ fl ],ilog[£t]} = | cot 

The nondegeneracy of the bracket follows immediately from {8 S , fit] = 2fj, t S st . Indeed, let / be 
a smooth function which is not constant. In that case, / depends nontrivially on at least one of 
the variables, say 6\, and hence Jjsfj- / on some open set. Therefore 

{/,Ai}=g^i,Ai} = 

on that open set. In other words, constants are the only functions which commute with every 
other function, and hence the Poisson bracket is nondegenerate. □ 

On this space, we are able to find the canonical coordinates for the Poisson structure defined 
above: 

Theorem 5. Define 

(4.4) q s = fi s JJ \z s - zj\, for all 1 < s < n. 
Then 

(4.5) {6 s ,q t } =2q t 5 st and {0 S , 9 t } = {q s , q t } = 0. 
In other words, 

(4-6) 0i, ... , 6 n , i log[qi], . . . , i \og[q n ] 

are canonical coordinates on the space of un-normalized measures fx. 

Proof. Having guessed the correct quantities which give the canonical coordinates, the proof is 
merely a calculation. As all the Zj's Poisson commute, the first bracket in (|4.5p follows immediately 
from P~2"j) . 

So we need to compute the bracket of the <7 s 's. For s^fwe have 
Us,qt} = {As, At} ' ]J \ z i ~ z& \ II ' Zfe _ z *' 

j^s kjtt 

+ a* n \ z 3 - z $\ • {As> n w - zt\\ 

j^s k& 

- a s n i zfe ~ z *i ■ ^ u n ~ z *\} 

fe#t j^s 

= ({fi s ,ft t } + T(s,t) - T(t,s)) ■ Y[ \zj - z s \ Y[ \z k - z t \, 
where, using ()4.2j) and the fact that s ^ t, we find that: 

rp( ,N ~ \ - {As, W ~ Z t \} 



2\z s -z t \ 2 

. _ ~ Z s Zt — Z s Zt 



({A s , z s }(z s - z t ) + {A 5 , z s }(z s - z t )) 
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By (|4.3p . the antisymmetry of T and using the fact the \zj \ = 1 for every 1 < j ' < n, we find that 

z s + z t , z sZt — Zs z t 



{p, s ,fit} + T(s, t) - T(t, s) = 2ipt, s p, t 
and hence 



%t \z s — z t 



2 



0, 



{q s ,qt} = o 

for any 1 < s, t < n. This proves the statement of the theorem. □ 

Once we have the result above, the canonical coordinates for the space of normalized measures, 
or, equivalently, of CMV matrices, follow from a simple observation: 

Corollary 4.2. With the notations from Theorem^ we get that 

(4.7) Ox,..., 0„_i, § log[ri, n ], . . . , | log[r„_ lin ] 

are canonical coordinates on the space of CMV matrices with fixed determinant, where 

f'j 



(4.8) r,, k = & [] 



zi - z k 



Proof. The observation that justifies our claim completely is 

1j 

r 3 ,k = — ■ 
qk 

Indeed, while each qj depends on the normalization through fij — c/ij, their ratios do not: 

Having observed this, the claim that (|4.7p are canonical coordinates follows by a moment's reflec- 
tion from l|4.5p . □ 

5. Compatible Poisson brackets 

In this section, we define a family of compatible (in the sense of Magri) Poisson brackets on the 
space of finite measures on the unit circle. Furthermore, the restrictions of all of these brackets to 
the manifold of probability measures represents a multi-Hamiltonian structure for the Ablowitz- 
Ladik equation (|1.2p described in the Introduction. 

Let ft be a smooth function on C which takes real values on S 1 . Define {•, - }h by specifying the 
bracket of the coordinates S and qi, 1 < s, t < n: 

(5.1) {0 s ,e t } h = {q s ,q t } h =O and {9 S , q t } h = 2h(e i9 »)q t S st , 

and then extend it in the canonical fashion: 

dfidf 2 dfxdf 2 



(5.2) {/i,/2}/, = X>(e 

s=l 



d9 s dq s dq s dQ s 



Proposition 5.1. Let h : S 1 — > K be a smooth, nonzero function. Then {-,-}h defined as in 
(15. 2[) ( or, equivalently, (|5.ip ) is a Poisson bracket, and, if h is not identically zero on any arc in 
S , then {-,-}h is nondegenerate. Furthermore, any two such brackets are compatible, in the sense 
that their sum is again a Poisson bracket. 

Remark 5.2. Note that for h identically equal to 1 we recover the extension of the Gelfand-Dikij 
bracket defined in Section 2J So the proposition claims that (|5.2p defines a family of Poisson 
brackets compatible with the GD-bracket. 

Proof. A simple calculation shows that {•, -}h obeys the Jacobi identity - note that it is sufficient 
to check it on the s 's and qt's. 

If h is not identically zero on any arc in S 1 , it follows immediately from the definition (|5.ip 
that the /i-bracket is nondegenerate: indeed, any smooth, nonconstant function / must depend 
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nontrivially on at least one of the variables, say 6\. Then the ^-bracket of / with the conjugate 
variable (in this case q\) will be nonzero: 

{/,<zi U - ^{hMh = 2h ( et9l ^§f 1 * °- 

So the only Casimirs are constant functions. 

Finally note that the newly-defined brackets are linear in h. In other words, for any h\ and fi2 
as above, the sum 

{"j + {") -}h 2 = {■) -}/n+/i 2 

is, by the previous argument, also a Poisson bracket. This is exactly the definition of compatibility. 

□ 

Going back to the fbj variables, direct calculations show that: 

Lemma 5.3. In the notations used above, the h-bracket can be written in the 9 and ft coordinates 
as: 

{z s ,jl t } h = 2iz s h(z s )fi t S st 

and, for s ^ t, 

z s + z t 



{/i s , in}h = ifrsfk (h(z s ) + h(z t )) 



z t 



The analog of the Carathcodory function for unnormalized measures is defined, unsurprisingly, 
by any of the following expressions 

/■ 1 n 
C + z A ~ lt \ _ sr^ ~ z i + z 



(5.3) F(z)= ^- d m = Y,^T^-= cF ^ 



S 1 C z pi Zj-Z 

where, as before, c = |/t| is the total weight of the finite measure jl, and F is the usual Caratheodory 
function associated to the probability measure /1 = fl/\jj\. Given the approach we take in this 
paper, it is natural to try to compute the ft-bracket of F at two distinct points z and w in the 
complex plane. The calculation that will give us these formulae is straightforward enough, with 
the only caveat that the resulting formula will involve not only F, but also the function 

F h ( Z ) = j2H^ 3 z f ±1 = jh{Q^dm- 

. =1 Zj Z J (, z 

As before, Fh=i(z) = F{z). Then one has 

{F(z),F(w)} h - i^^{F(z) - F(w)) (F h (z) - F h (w)) 

- iF(0)(F h (z) - F h (w)) - iF h (0)(F(z) - F(w)) 



.A (zj + z k )(h(zj) + h{z k ))(zjZ k + zw) . _ 
- t(W {zj - z){zk - z){zj - w){zk - w) W 



While this formula in fairly involved and not very pretty, it simplifies greatly when restricted to 
the manifold of probability measures. In order to find this restriction, we need to take the reverse 
road to that in the previous section, and hence compute the following bracket: 

r £./ u . f (zj + z k )(h(zj) + h(z k )) . . 
{c,F(*)} fc = ,*£ {z ._ z){zk _ z) W 

= i[F h (z)F(z) - F h (0)F(0)] 

{c,F{z)} h = ic[F h {z)F{z)-F h (0)} 
(where Fh{z) — cFh{z), and recall F(0) = 1). Finally, we obtain the following: 
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Theorem 6. The restrictions of the h-brackets to the manifold of probability measures supported 
at n points on the unit circle are given by 



(5.4) 



{F(z), F(w)} h = i(F h (z) - F h (w)) ^ z {F(z) - F(w)) + F(z)F(w) - 1 



which defines, for h smooth and real valued, a family of compatible Poisson brackets, that forms a 
multi-Hamiltonian structure for the defocusing Ablowitz-Ladik bracket. 
In particular, this implies 

(5.5) {z s ,z t } h = 0, {z s ,fi t } h = 2iz s h{z s )nt{8 s t - Ms) 

and, for s ^ t, 

r i . [(h(z s ) + h(z t ))(z s + z t ) (h(z t ) + h{z k )){z t + z k ) 

{H s ,Ht}h = t^sfJ-t > Mfc 1 

^-^ L z s — zt 

(5.6) 

(h(z k ) + h(z s ))(z k + z s 



zt - z k 



Z k Z s 

Proof. The equation (|5.4|) follows directly from the previous formulae and from 

{F(z),F(w)} h = c 2 {F(z), F(w)} h + cF(z){c, F(w)} h + cF(w){F(z), c} h . 

The formulae for the brackets of the z s 's and /xt's follow from (]5 .4[) by using the residue theorem, 
exactly as in the proof of Theorem 2J □ 

Finally, we close this section by identifying Hamiltonians for Ablowitz-Ladik flows in the h- 
brackets. 

Proposition 5.4. Let g be a polynomial, and consider the Hamiltonian on M n (C) defined by 



4>{X) = ImTr(<?(X)). Then the evolution of the spectral 



measure 



associated to a CMV matrix C under this Hamiltonian in the h-bracket is given by 
(zj ={<j,, Zj } h = 

\Aj = {<t>,V>j}h = Vj G{zj)h{ Zj ) - J2?=i G(zi)h(zi)m , 
where G(z) = 2Re(zg'(z)). Equivalently, 

e G(z)h(z)t^ t = o) 



(5.8) 



\ e G{z)h(z)t^ t = 0)| 



Proof. Note that we can rewrite (j> = Im£fc = i g(z k ). Then the first formula in (|5.7|) follows 
immediately from (|5.5[) . while 



= {^Vj}h =^Qg-{0k,fJ,j}h 



k=l 



But 



and hence 



^=Im(ie i V(e^)) = |G(e^), 



MfeJ 



k=l 



n 

^[G(e^)h(e^)-J2G(e^)h(e* e <)m 



1=1 



as in (|5.7[) . Finally, (|5.8p follows from (|5.T|) by integration. 



□ 



Recall that the Gelfand-Dikij bracket, which is the Poisson bracket associated to the Ablowitz- 
Ladik equation, corresponds to h = 1, while the Hamiltonians which generate the flows in the 
AL hierarchy are exactly of the form </>(C) = ImTr(g(C)) for a polynomial g. So we find that the 
following holds: 
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Corollary 5.5. If h is a trigonometric polynomial, then the flow generated in the h-bracket by 
the Hamiltonian 4>{C) = ImTi(g(C)), where g is a polynomial, is one of Ablowitz-Ladik flows. 

Proof. Since both h and G are trigonometric polynomials, their product will have the form 

d 

h(z)G(z) = Cjjy 7 ', where c^j—Cj. 

J = -d 

Then for z E S 1 

d 

h{z)G{z) = 2Re(J2c 3 z j ) + c = 2 Re(zg'(z)) + c Q = G(z) + c , 

3=1 

where we can choose g(z) = Y]j—i (unique up to an additive constant). If we set 0(C) = 
ImTr(<7(C)), then we see that 

n 

{<P, t*j}h ~ {4>, Mj}i = fij [ c o ~ ^2 Co ^] = °' 

k=l 

since [i is a probability measure. This is exactly what we claimed. □ 



We close this section by noting an immediate consequence of Corollary 15.51 Consider two 
trigonometric polynomials, hj, j = 1,2, and two polynomials gj from which we construct Hamil- 
tonians as before: <j>j{C) = lmTr(gj(C)). Then 4>i generates the same flow in {•, as </>2 does in 
{•,•}/* iff 

(5.9) hiGt - h 2 G 2 = const, e R, 

where, as before, Gj(z) = 2Re(zt^(z)) for j = 1,2. 

6. The connection to Schur flows and the Toda lattice 

Let us now consider the case where the measure dfi is symmetric with respect to complex 
conjugation, or what is equivalent, where all the Verblunsky parameters are real. In this case, 
there are an even number of eigenvalues, z\, z n , n = 2N, with the extra symmetry 

(6.1) Zj+N = Zj, [ij+N = Lij for l<j<N. 

For simplicity of the notation, we further assume that z\, z/v are the eigenvalues on the upper 
half of the unit circle. It is a famous observation of Szego (see [Szel §11.5]) that the polynomials 
orthogonal with respect to this measure are intimately related to the polynomials orthogonal with 
respect to the measure dv on [—2, 2] defined by 

(6.2) f f(z + z- 1 )d f i(z)= f f{x)dv{x). 

JS 1 J-2 

The recurrence coefficients for these systems of orthogonal polynomials are related by the Geron- 
imus relations: 



(6.3) 



jbk+l = (1 - U 2 k-l)et2k ~ (1 + Ci 2 k-l)ci2k-2 
1 1/2 

[ak+i = {(1 - a 2 fc-i)(l - al k )(l + a 2 k+i)} 

It is an easy observation (see, for example, |Nen2p that the second flow in the Ablowitz-Ladik 
hierarchy, which is generated by Im TrC, will preserve the property of all the Verblunsky coefficients 
being in (—1,1). Hence it makes sense to ask what is the flow it induces via the Geronimus 
relations to the a's and 6's. A direct calculation shows that the answer is exactly the Toda flow! 
It is immediately clear from (|1.5I) that the submanifold of real Verblunsky coefficients, which we 
will denote by M, is in fact stable under any of the flows generated in the usual AL (or GD) 
bracket by the Hamiltonians 1th Kk = i ImTr(C fc ), k > 1. We will call these the Schur flows ( cf. 
FavGckl] IGol| . But neither the Hamiltonian in question, nor the usual Ablowitz-Ladik Poisson 
bracket have meaningful restrictions to this submanifold, nor is it possible to find the image of the 
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flows generated by Imifj. under the Geronimus relations by straightforward calculations. In this 
section we investigate the newly defined /i-brackets from these points of view. 

Proposition 6.1. The h-bracket {v}/i has a restriction to the submanifold M of probability 
measures on the unit circle which are symmetric with respect to complex conjugation iff h(z) — 
-h{z) for z G S 1 . 

Proof. To prove this statement, consider a set of functions on the manifold of measures supported 
at n = 2N points defined by 

(6.4) X s — z s + z s and v s = 2\x s - 

In fact, these functions with 1 < s < N form a set of coordinates on the submanifold of probability 
measures on the unit circle which are symmetric with respect to complex conjugation, and 

(6.5) A s = z s + z s+N , v s = /i s + fi s+N 

for 1 < s < N. Furthermore, the measure dv on [—2, 2] defined above is exactly dv = X^fLi v s^\, ■ 



Direct calculation using (|5.5[) shows that 

{A s , At}/j = 0, for every 1 < s, t < N. 

The other two types of brackets are more complicated. Decompose a general function h as h = 
h+ + h-, where h + (z) — h+(z) and h^(z) — —h-(z) for z £ S 1 . This is equivalent to setting 
2h + {z) = h(z) + h(z) and 2h-(z) = h{z) - h(z). Then consider 1 < s,t < N, and use |63|) . (|6~Tj) , 
and (|5.5p with s,s + N, t, and t + N, respectively: 



{K,vt}h = 2iz s h(z s )^ t (S st - Us) + 2iz s h(z s )nt+N(-Hs) 

+ 2iz s+N h(z s+N )ii t (-^, s+N ) + 2iz s+N h{z s+N )^ t+N {8 st - ^ s +n) 
= h + (z s ) [2i(z s + z s )[i t (6 st - Us) - 2i(z s + z s )^ s ^t] 

+ h-(z s )[2i(z s - z s )pL t {5 s t - fJ-s) - 2i(z s - z s )pL s pL t ] 
= 2i(z s + z s )h + (z s )fi t (S st - 2fj, s ) 
+ 2i(z s - z s )h-{z s )^t{8st - 2pL s ) 
The last expression is real-valued iff h+ = 0, or, equivalently, h{z) — —h(z). In this case, we get 

{X s , v t } h = 2i(z s - z s )h^(z s )nt(5 st - 2^i s ), 

where the right-hand side is real for z s £ S 1 , and invariant under the mapping taking a probability 
measure to its complex conjugate. 

Finally, we need to deal with the ^.-bracket of the ^'s. Proceeding as above, we get: 

{ v si v t}h = {Ms> Mt}/i + {^s+N, Mt+Jv}/i + i^s+N, f-t}h + {Us, Mt+Jv}ft, 

Group the first two terms on the right-hand side to get 

-(h(z s ) + h(z t ))(z s + z t ) (h(z t ) + h(z k ))(z t + z k ) (h(z k ) + h(z s ))(z k + z s )- 



k=£s,t 



+i^s+N^t+N ^ Mi 
l^s+N,t+N 



Zt Z t - Z k Z k 

(h(z s+N ) + h(z t+ N))(z s +N + z t+N ) (h(z t+N ) + h(zi))(z t+ N + zi) 

Z s +N — Z t +N Z t +N — Zi 

(h(zi) + h(z s+N ))(zi + z s+N ) 



k^s,t 



Zl — Z s +N 

(h(z 8 ) + h(z t )) (z s + z t ) (h(z s+N ) + h(z t+N )) (z s+N + z t+N ) 



Z s — Zt Z s +N — Zt+N 

cyclic permutations 
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To obtain this last identity, set I — k + N in the second sum on the left-hand side, and use formula 
(15. 6p as well as the symmetry conditions (|6.ip . Note that we think of the indices as periodic, with 
period 2N. The main observation at this point is that, for any two indices 1 < Z, r < N, we have 

(h(zi) + h(z r )) (zi + z r ) (h(zi +N ) + h{z r+N )) (z l+N + z r+N ) 



Z\ — Z r Zi + n — Z r+ N 

= (^ + (z i ) + / l+ (z r ))[fi±i!l + fi±ilj + (h-(zt) + h-(Zr)) 
But it is a simple observation that, for z, w G S , 

Z + ill Z + w 

z — w z — w ' 

and hence the expression above equals 

Zl + Z r 



zi + z r zi 



Zl - Z r Zl 



2{h-{zi) +h-(z r ))- 



zi - z r 

In particular, this means that we can work our way backwards to the original expressions for the 
four /i-brackets, only with h replaced by h-. Note that, for any z ^ w G S , ijz^ € K and 

(h-(z) + h-(w))(z + w) (h-(z) + h_(w))(z + w) 
z — w z — w 

But this immediately implies that, even though the bracket {v s ,^t}h is still the sum of four 
complicated formulae, it is real-valued and invariant under the mapping taking a probability 
measure to its complex conjugate, which completes the proof. □ 

Now consider a function h obeying 

h{e- ie ) = -h{e ie ), 

and restrict the bracket {•, -}h to the subspace of real Verblunsky coefficients. We want to write 
this restriction as a combination of the (compatible) Poisson brackets for the Toda lattice found 
in FayGck2 . To avoid confusion, we will go back to the notation from Section [3] Hence we will 
start denoting the /i-brackets by {•, , the superscript denoting the fact that this is a Poisson 
bracket compatible with the Gelfand-Dikij bracket. By contrast, we will later on be interested 
in some of the Poisson brackets compatible with the Lie-Poisson bracket, which were originally 
introduced in |FayGek2] , and which we will denote here by {•, for some function H. 

In order to achieve this, we must first relate the Caratheodory function of a measure dpi on the 
circle which is invariant under complex conjugation to the m-function of the associated measure 
dv on [—2, 2]. A simple calculation (see, for example, [Simlp show that, for z G C \ {0}, we have 

(6.6) F{z) = -F^) = {z-\)-m{z + \), 
and 

(6.7) F h (z) = F h (l)=im H (z + l), 
where 

(6.8) H(2cos8) = 2sm{6)h(e l9 ), 
h obeys h{e~ 10 ) — —h(e lS ) as above, and ran is defined by 

f 2 1 

(6.9) m H (X) = J H{t)j-^d v {t). 

Note that, for z G S 1 , H is defined by 

H(z + z) = (z- z)h(z), 

which is exactly the type of combination that has already appeared in the proof of Proposition ^. II 
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We wish to find the ^.-bracket of the m function and compare it to the ones in |FayGek2] . To do 
so, we insert (|6.6[) into (|5.4p . but we must not lose sight of the symmetry inherent to the situation. 
In this case, we use the fact that 

(6.10) 4(z -!)(«,- i){m(z + ±),m(™ + = {F(z) F(±),F(w) 

Using (|6.6p . we get that 

{f W - F(l),f (») - = {FW,f H }W . {f (W )}W 



(i) 



where 



B(z, w) 



= (F h (z) - F h (w)) ■ (B(z,w) - 4F{z)F{w)) 
w 4- - 

(F(z) - F(w)) + z -{-F{z) - F(w)) 



w — z 

i + Z 

+ f (F(z) + F(w)) - f f(- F ( z ) + F ( w )) 



w ■ 

z 

1 i 1 



= 4z 



— z 

l-w 2 



{z — w)(l — zw) 



F(z) - Aw 



1-z 2 



(z — w)(l — zw) 



F{w) 



, . , , . . m(z + - ) — rn(w + — ) 

= 4(*-;)(<«-£)- ' " 



Putting it all together and replacing the combination z + ^ by a general A £ C \ spec(<iz/), we 
get that 



{m(A),m(0}£° = (m H (\)-m H (0) 



m(A) — to (£) 



m(A)m(^) 



Let us recall that, for fc > 0, Faybusovich and Gekhtman defined a Poisson bracket on the 
space of measures supported on [—2,2]. This fc-bracket can be written down in terms of various 
coordinates, but here we concentrate on its expression in terms of the associated m-functions: 

'm(X) — 77i(£) 



(6.11) {m(A),m(£)} fe = ((A fc m(A))_ - (£ fc m(0)-)) 



— m(A)rn(£) 



0, 



where, for a meromorphic function r : C - 

AT 

MA))+=5> A ' and (r(A))_=r(A)-(r(A)) + = £ qA 



with r(A) = X);=-oo C 'A Z ' iV < oo, we set 

-l 



z=o 



i=-c 



Note that for any k > we have 



t fc - A 
t-X ) 

x k 



k — 1 



= (^Va*-*- 1 ) =0 



t — A 



and hence 

(t- a)- = it- xJ- 

So by integrating, we get that 

(H(X)m(X))_ =m H (X), 

where H is a polynomial and m# is defined by (|6.9[) . In other words, we have proved the following 

Proposition 6.2. If h is a smooth, real-valued function on the unit circle S 1 such that h(z) — 

—h{z), then the restriction of the bracket {•, to the submanifold of probability measures in- 
variant under complex conjugation is given, for A ^ £ C \ K, by 



(6.12) 



{m(A),m(0}i 1) = (m H (A)-m ff (0) 



77T,(A) — 77l(£) 



m(A)m(£) 



^We use here the notation of [FayGek2] for the compatible brackets. 
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where H , m and mjj are defined as above. 

The right-hand side defines a Poisson bracket {-,-}^P on ^ e m o>nifold of probability measures 
supported at N points on [—2,2] which is compatible with the Toda lattice (i.e. restriction of the 
Lie-Poisson) bracket. Furthermore, if h is a trigonometric polynomial, then H is a polynomial 
and {•, is a linear combination of the compatible brackets (|6.11[) of |FayGek2] . 

While we could try to write down the general formula for these brackets, we will limit ourselves 
to investigating the simplest case, which is 

h{e ie ) = 2sra(0). 

Then we get that H(2cos(8)) = 4sin(6') 2 , or, equivalently, 

H(t) = A-t 2 , for te[-2,2]. 

Therefore we find that the restriction of the Poisson bracket {•, ^^sin ^° * ne s P ace °f rea l Verblunsky 
coefficients coincides, under the Geronimus relations, to the Poisson bracket {•, -yg — 4{-, -} — 
{;-h- 

We close this paper by identifying the Hamiltonians defining certain Ablowitz-Ladik Hows with 
Toda hierarchy Hamiltonians. Recall that the usual AL bracket is {•, ■}f ll for hi = 1, and let 
h 2 (z) = — h 2 (z) for z G S 1 , as in Proposition [ITU Consider a polynomial with real coefficients, gi, 
and note that the Hamiltonian <p\{C) — ImTrgi(C) is a linear combination, with real coefficients, of 
the flows ImKk and hence it generates a flow under which M is stable. If g 2 is another polynomial 
and 4> 2 (C) = ImTrg 2 (C), then from (j5l)]) we get that if 

(6.13) GiO) - h 2 {z)G 2 {z) = const, G R, 
then the two flows coincide: 

{01: -} G D = ih, -}h 2 ■ 

Here, as in Section^ Gj(z) = 2Re(z^(z)). Let us make a few remarks on (16. 13[) : 

• Since gi(z) — X)j=o c i z "' w ith Cj £ R, we get G\{z) — 2j=i J c j( z "' + ^) an d so for z G S 1 , 
G\(z) = G\(z). Since we imposed ^2(2) = —h2(z), we must have that G2(z) = —G2(z), or, if we 
work as in the proof of Corollarv l5.5l the coefficients of 52 are purely imaginary, up to an additive 
constant which we will ignore since it does not influence the flow. So the first observation is that 

(6.14) 4> 2 (C) = ReTr(-i 92 (C)) 

is a linear combination with real coefficients of the Hamiltonians Reiffc, and hence has a nontrivial 
restriction to the submanifold M. 

• A straightforward count of the parameters in (|6.I3p shows that it is not true that given gi, we 
can always find h 2 and G 2 with the required properties and which satisfy (|6.I3[) . Indeed, without 
loss of generality we may assume that g\, h 2 and G 2 are monic. Then the right-hand side of (|6.f 3p 
is determined by deg(/i 2 ) — 1 + deg(G*2) — I = dcg(Gi) — 2 parameters, while the left-hand side 
imposes deg(Gi) — 1 conditions. In other words, we need an extra degree of freedom. Without 
going into too many details, let us mention that one way to deal with this problem is to allow 
the Schur flow in question to be modified by a constant multiple of the first Schur flow: for any 
gi monic of degree at least 2 and with real coefficients, there exist a real constant c, and monic 
h 2 and g 2 as above so that G\(z) = 2Ke(zg[(z) — cz), h 2 and G2 obey (|6 . 13[) . In particular, the 
restriction to the submanifold M of a Schur flow whose Gi obeys f|6 . 1 3[) for some appropriate h 2 
and G2 is a Hamiltonian flow in the Poisson bracket {•, •} \mU 

• Finally, we must understand the image through the Geronimus relations of the restriction to 
M of Hamiltonians given by Re(g(C)), with g polynomial with real coefficients. But this can be 
obtained immediately from Proposition B.3 of Killip and Nenciu [KilNenl] . which shows that, if 
ctj G (—1,1) for all j, then C + C T is unitarily equivalent to a direct sum of two Jacobi matrices, J 
and J, where entries of J are related to the Verblunsky coefficients defining C by the Geronimus 
relations. Furthermore, the spectral measure for J w.r.t. the vector e± = [10 - ■ - ] T G R™ is dv, 



6 Since h 2 (z) = 



—h,2(z), we know from Proposition ^, ll that M is a Poisson submanifold in this Poisson structure. 
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while the spectral measure of J w.r.t. e\ is dv(x) — 2 (i-a'^)(i- ai ) — x 2 )dv(x). In particular, 

we see that the two measures v and v have the same support, or, equivalently, J and J have the 
same eigenvalues. For example, 

ReTrC = ±Tr(C + C T ) = ±Tr(J © J) = Tr(J). 

A slightly more careful analysis will show that for any k > 1 there exists a monic polynomial g^, 
with real coefficients, so that 

ReTr(C fc ) \ M = Tr( 3fe (J)). 
But the Hamiltonians on the right-hand side are exactly Toda hierarchy Hamiltonians. 
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